X-ray Observations of the Young Pulsar J 1357— 6429 and Its 

Pulsar Wind Nebula 
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ABSTRACT 



We observed the young pulsar J1357-6429 with the Chandra and XMM-Newton observatories. 
The pulsar spectrum fits well a combination of absorbed power-law model (r = 1.7 ± 0.6) and 
blackbody model [kT — 140^4$ eV, R ~ 2 km at the distance of 2.5 kpc). Strong pulsations 
with pulsed fraction of 42% ± 5%, apparently associated with the thermal component, were 
detected in 0.3-1.1 keV. Surprisingly, pulsed fraction at higher energies, 1.1-10 keV, appears to 
be smaller, 23% ± 4%. The small emitting area of the thermal component either corresponds to 
a hotter fraction of the neutron star (NS) surface or indicates inapplicability of the simplistic 
blackbody description. The X-ray images also reveal a pulsar-wind nebula (PWN) with complex, 
asymmetric morphology comprised of a brighter, compact PWN surrounded by the fainter, much 
more extended PWN whose spectral slopes are T = 1.3 ± 0.3 and T — 1.7 ± 0.2, respectively. 
The extended PWN with the observed flux of ~ 7.5 x 10~ 13 erg s _1 cm~ 2 is a factor of 10 more 
luminous then the compact PWN. The pulsar and its PWN are located close to the center of 
the extended TeV source HESS J1356-645, which strongly suggests that the VHE emission is 
powered by electrons injected by the pulsar long ago. The X-ray to TeV flux ratio, ~ 0.1, is 
similar to those of other relic PWNe. We found no other viable candidates to power the TeV 
source. A region of diffuse radio emission, offset from the pulsar toward the center of the TeV 
source, could be synchrotron emission from the same relic PWN rather than from the supernova 
remnant. 

Subject headings: pulsars: individual (PSR J1357-6429) — X-rays: individual (CXOU J135605. 9-642909, 
1RXS J135605. 5-642902) — ISM: individual (HESS J1356-645, G309.8-2.6) 



1. Introduction 

Thanks to their rich observational manifesta- 
tions, young energetic pulsars are among the most 
attractive targets for the Chandra and XMM- 
Newton observatories. These manifestations in- 
clude thermal emission from the neutron star (NS) 
surface, non-thermal emission from the pulsar 
magnetosphere, and a pulsar-wind nebula (PWN) 
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produced by the interaction of the pulsar wind 
with the ambient medium. The magnetospheric 
and thermal emission can differ in their rela- 
tive strengths and exhibit distinct spectra and 
pulse profiles. The magnetospheric emission usu- 
ally has power-law (PL) spectrum and is strongly 
pulsed. Detailed studies of several bright young 
pulsars indicate a nonuniform NS surface tem- 
perature distribution with cool (Tc < 100 cV) 
and hot (Th ~ 0.1 — 0.3 keV) thermal compo- 
nents, which are believed to be emitted from the 
bulk of the NS surface and polar caps, respec- 
tively (e.g. JPavlov et al.1 l2001i I Zavlin fc Pavlov 



2004 : Ide Luca et al.ll200a iKargaltsev et al.ll2005T 
Since the sample of young pulsars with high qual- 
ity X-ray spectra is still very small, further obser- 
vations of young pulsars must be carried out to 



learn about the early epoch of NS cooling and test 
magnetospheric emission models. 

In addition to X-ray emission from the NS sur- 
face and magnetosphere, a significant amounts of 
non-thermal emission is radiated from the PWN. 
PWNe display various shapes and structures, in- 
cluding (but not limit ed to) jets and tori, o r 
cometary shaped tails (jGaensler fc Sland 120061 ) . 
Observations with high angular resolution are par- 
ticularly useful since they resolve the extended 
PWN emission, seen downstream of the termina- 
tion shock, from the magnetospheric non-thermal 
emission. Comparing properties of these two com- 
ponents, one can learn about the physical pro- 
cesses responsible for an efficient but yet un- 
known particle acceleration mechanism energizing 
PWNe. 

Although more than ^70 PWNe have been 
observ ed with Chandra and XMM-New ton to 
date (JKargaltsev fc Pavlovl 120081 |2010L here- 
after KP08 and KP10), the existing models 



( Kennel fc Coronitill984i Komissarov fc Lvubarskv 



2004; Ivan der Swaluwl l2005t iBucciantini et al 
20051 ) still cannot explain some of the observed 
complex morphologies. Pulsars of similar ages 
and spin-down luminosities can have PWNe with 
quite different properties. As of now, it is not clear 
whether these differences can be solely attributed 
to different properties of the surrounding interstel- 
lar medium (ISM) or different intrinsic properties 
of pulsars also play an important role. A larger 
sample of well-resolved X-ray PWNe associated 
with pulsars having different intrinsic properties 
and residing in different environments is needed 
to further advance our understanding of pulsar 
winds. 

The young pulsar PSR J1357-6429 (hereafter 
J1357; P = 166 ms, E = 3.1 X 10 36 erg s"\ 
surface magnetic field B = 7.8 x 10 12 G) was 
discove red in the Pa r kes M ultibeam Pulsar Sur- 
vey by ICamilo et al.1 (J200J). The spin-down age, 
t = P/(2P) = 7.3 kyr, indicates that J1357 is one 
of the youngest pulsars known. The distance of 2.5 
kpc, estimated from the pulsar's dispersion mea- 
sure (DM = 127.2 cm" 3 pc) and the Galactic free 
electron density model (|Cordes fc Lazioll2002l ). im- 
plies a spin-down flux E/(And 2 ) = 1.2 x 10~ 9 erg 

— 1 —2 

s cm . 

The pulsar is located within the T eV source 
HESS J1356-645 (hereafter HESS J1356: lRenaud et al 



20081 : lAbramowski et al.l 120111) . The extended 
TeV emission is detected up to 12' from the 
peak of the TeV surface brightness, located at 
R.A.= 13 h 56 m , deel- -64°30'. Given the 
relatively small angular separation, ~ 7', be- 
tween J1357 and the center of the HESS source, 
and the TeV luminosity, LtcV ~ 6x 10 33 ergs 
s ~ 0.002S (at the 2.5 kpc distance), it seems 
plausible that the HESS source could be powered 
by this pulsar. The TeV spectrum of HESS J1356 
fits PL with the photon index TtoV ~ 2.2, without 
a cutoff at high energies. A possible detection of 
J1357 with AGILE at lower 7-r ay energies (100 
MeV - 10 GeV) was reported bv IPellizzoni et al 



(|2XT09jjl|. No optical counterpart t o the pulsar or 



its P WN has been detected yet (jMignani et al 



20111) 



At radio frequencies, iDuncan et al.l (|1997l ) re- 
ported the extended emission from the SNR can- 
didate G309.8-2.6 , which is spatially coincident 
with HESS J1356. The extent of the radio emis- 
sion appears to be smaller than that of the TeV 
source. J1357 has been previously observed with 
XMM-Newton (2005 August 5; 15 ks; PI F. 
Camilo) and Chandra HRC-S (2005 November 
18 and 19; 16 and 17 ks; PI M. Mendez, ob- 
server L. Kuipe r). The ana lysis of those data 
was re ported by IZavlinl (|2007j ) and lEsposito et al. 
(J2007l l 2 l. These authors found that the pulsar 
spectrum likely c onsists o f ther mal and nonther- 
mal components. IZavlinl (J2007) found pulsations 
with the radio pul sar period and pulsed fraction of 
63% ± 15%, while lEsposito et al'l (J2007I ) reported 
only an upper l i mit of 30% on modulation am- 



plitude. IZavlinl (J2007I ) also reported a tail-like 



extended emission in the HRC-S image, with a 
lumin osity of 2.5 x 10 3 1 erg s -1 (for d = 2.5 kpc), 



while lEsposito et al.l (|2007l) concluded that the 



1 After our paper was submitted, iLemoine-Goumard et aD 
(2011) reported results of observations of J1357 with Fermi 
LAT. They detected 7-ray pulsations and measured the 
photon index Tg c v ~ 1.5, found an exponential cut-off at 
~ 0.8 GeV, and estimated the luminosity I/GeV ~2x 10 34 
ergs" 1 , for E > 0.1 GeV. 

2 It should be noted that first report on the Chan- 
dra results was published in "Chandra News" 
(March 2007, Issue number 14, pages 12-13 (see 

http : //exc . harvard . edu/newsletters/news_14/newsletterl4 . html ) 
by R. Kraft and A. Kenter. In particular, evidence for the 
presence of extended emission (PWN) and the detection of 
pulsed signal (pulsed fraction 40% ± 12%) were reported 
in that publication. 



surface brightness distribution around the pulsar 
"is consistent with that from a point source" and 
estimated a 2-10 keV upper limit of 3 x 10 31 erg 
s _1 on the PWN luminosity. To resolve these 
controversies, reliably disentangle various emis- 
sion components, and study the relation between 
the pulsar and HESS J1356, we have carried out 
deeper observations with the Chandra ACIS and 
XMM-Newton EPIC detectors. 

The details of the observations and data analy- 
sis are presented in Section 2. We discuss possible 
interpretations of the PWN morphology, describe 
inferences from the pulsar spectrum, and specu- 
late on the nature of HESS J1356 and its relation 
to the other sources in the field in Section 3. 

2. Observations and Data Analysis 

We observed J1357 for 60 ks with the ACIS 
detector on board Chandra on 2009 October 8 
(ObsID 10880). The target was imaged on the 
13 chip, ss 30" from the aim point. The other 
chips operated during the observation were 10, II, 
12, S2, and S3. The observation was carried out 
in timed exposure mode and telemetered in Very 
Faint format. The detector was operated in full 
frame mode, which provides time resolution 3.24 
s. The useful effective exposure time (live time) is 
59.22 ks. There were no significant particle back- 
ground flares during the observation. The data 
were reduced and analyzed with the Chandra In- 
teractive Analysis Observations (CIAO) package 
(ver.4.2), with CALDB 4.1.4. We used Chandra 
Ray Tracer (ChaRTJl and MARX softwareQ for 
image analysis and XSPEC (ver.12.5.0) for spec- 
tral analysis. 

J1357 was also observed with the XMM- 
Newton EPIC MOS1, MOS2 and PN detectors on 
2009 August 14 and 15 (ObsID 0603280101). The 
two MOS detectors were operated in Full Window 
mode with a time resolution 2.6 s, while the PN de- 
tector was in Small Window mode, which provides 
6 ms time resolution (at the expense of reduced 
efficiency, « 71%). The observation was processed 
with the XMM-Newton Science Analysis Software 
(SAS; ver.9.0.0). The total live times of observa- 
tions were 78.1 ks and 55.5 ks for EPIC MOS and 




3 See http://cxc.harvard.edu/chart/threads/index.html 
4 See http://space.mit.edu/CXC/MARX/. 



Fig. 1. — Chandra images of J1357 and its PWN 
in the 0.3-8.0 keV band. The exposure-corrected 
larger image is binned to a pixel size of 0'.'98 
and then smoothed with the Gaussian kernel of 
about 3" radius using the aconvolve script. The 
dashed lines show the chip boundaries. PWN II 
and PWN 12 are the extraction regions (see text). 
The zoomed unsmoothed image in the inset shows 
the extraction regions for the image and spectral 
analysis. 



PN, respectively. However, the data were con- 
taminated by strong flares, which exceeded the 
quiescent count rate (0.1 counts/s for MOS and 
0.2 counts/s for PN) by a factor of up to 8. After 
excluding contaminated data, the useful scientific 
exposure times are 62.2 and 44.6 ks for MOS and 
PN, respectively. The medium filter was used for 
all the EPIC detectors. The data were filtered to 
allow only standard event grades (patterns < 12 
for MOS and < 4 for PN). 

2.1. Images 

2.1.1. Chandra 

To produce Chandra images of the pulsar and 
its vicinity at subpixel resolution, we removed 
the pipeline pixel randomization and applied 
the subpixelization tool which improves quality 
of the images using grades of split-pixel events 
(|Mori et al.lboOll ; iTsunemi et al-lboOlh . Figured] 
shows the Chandra image of the region around 
J1357 and the extraction regions on the ACIS- 
13 chip used for the image and spectral analy- 
sis in the 0.3-8 keV energy band. In the im- 



age we see a bright pointlike source surrounded 
by diffuse emission. The position of the center 
of the brightest pixel is R.A. = 13 h 57 m 02!496, 
dccl.=-64°29'30"06 (J2000.0). The radio pulsar 
positionisR.A.=13| 1 57 m 02 i ;43(2),decl.=-64 o 29 / 30"2(l) 
(jCamiloet alll2004l . The difference of 0"45 be- 
tween the X-ray and radio positions is smaller 
than the error in absolute Chandra astrometry 
(0'.'6 at the 90% confidence levels). We measure 
the background on the 13 chip in the 50" radius 
circle centered at about 1 '5 southeast from the 
source. The bright pointlike source contains 371 
counts in the r = 1" aperture, including 12% 
contribution from the background and the diffuse 
emission around the source. 

The diffuse emission in the immediate vicinity 
of the pulsar (region PWN1; shown in the inset 
in Figure [1}, excluding the region of r = 1'.'5 cir- 
cle centered on the pulsar, has 282 total counts in 
91 arcsec 2 , including 4.6% background contribu- 
tion. The background-subtracted surface bright- 
ness within the PWN1 region is ~ 2.97 ± 0.19 
counts arcsec -2 . We also define PWN2 region (see 
the same inset in Figured]), which includes the dif- 
fuse feature extending northeast of the pulsar and 
turning northward at ~ 23" from the pulsar. This 
region contains 74 counts in 108 arcsec 2 , including 
21% background contribution. This corresponds 
to the background-subtracted surface brightness 
of ~ 0.54 ± 0.08 counts arcsec -2 . 

The larger exposure-corrected, binned and 
smoothed image (Figure []} shows an even fainter 
diffuse emission extending ~ 2. '5 westward and 
over 1.'5 northward. The PWN region on II 
chip (PWN II; wl921 arcsec 2 ) includes 295 to- 
tal counts, of which about 70% comes from the 
background, estimated from the 150" radius circle 
at about 5' north-northeast of the source region. 
The total number of counts in the PWN region 
on 12 chip (PWN 12; «6128 arcsec 2 ) is 1036, in- 
cluding 72.5% contribution from the background, 
estimated from the 100" radius circle located at 
just south of the source region. The background- 
subtracted PWN surface brightness within the 
source regions on II and 12 chips is 0.046 ± 0.009 
and 0.047 ± 0.006 counts arcsec -2 , respectively. 

The radial profile of the emission centered on 
the brightest pixel position is plotted in Figure [2] 




5 See |http://cxc.harvard.edu/cal/ASPECT/celmon/| 
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Fig. 2. — Radial profiles of the observed emis- 
sion (solid histogram with bin size 0"25) around 
the pulsar and the simulated PSF (dashed line) in 
the 0.3-8.0 keV energy band. The PWN emission 
dominates the pulsar PSF at r > 0"75 from the 
pulsar. 

with the simulated point spread function (PSF) 
shown for comparison. To simulate the PSF, we 
used the ChaRT and MARX packages following 
the standard procedure. To reduce the statisti- 
cal errors of the simulation, the normalization of 
the input spectrum was increased by a factor of 
100 while running ChaRT. The simulated images 
were then rescaled back for direct comparison with 
the data. We have also tried several values of the 
MARX Dither Blur parameter and found an op- 
timal value of 0"25, which minimizes the differ- 
ence between the simulated PSF and the observed 
count distribution in the immediate vicinity of the 
pulsar position. The histogram in Figure [5] shows 
that the PWN starts dominating the unresolved 
pulsar emission at > 0"75 from the pulsar. 

2.1.2. XMM- Newton 

Figure|3]shows the combined EPIC MOS1+MOS2 
image (in 0.3-10 keV) and several regions that we 
use in our analysis. Within the r — 20" circle cen- 
tered on the pulsar we find 1,900 MOS1+MOS2 
counts, of which the background (taken from the 
r = 100" circle east of the pulsar) contributes 
about 15%. The MOS images show the large-scale 
PWN morphology reminiscent of that seen in the 
Chandra image (cf. Figure [1} . Indeed, the con- 
tours from the MOS1+2 PWN image match the 





Fig. 3.— XMM-Newton EPIC M0S1 + MOS2 
image of J1357 and its PWN in the 0.3-10 keV 
energy band. The image is binned to a pixel size of 
3" and then smoothed with a Gaussian kernel of 9" 
radius. Extraction regions (r = 20" circle for the 
pulsar region, 11210 arcsec 2 polygon for the PWN 
region, and r = 100" circle for the background 
region) for the spectral analysis are also shown. 

Chandra PWN image quite well (see Figure 0J. 
Since the source is near the chip edge in the PN 
image, only part of the extended emission can be 
seen, while the other part is out of the held of view 
(see Figure [5]) . We choose anr = 20" circle cen- 
tered on the pulsar in the PN image for the source 
region, and a 60" radius circle, centered at about 
2. '5 south-southwest from the source for the back- 
ground region. The total number of counts in the 
source region is 2,403, including 19% background 
contribution. 

To choose an optimal PWN region in the XMM- 
Newton MOS and PN images, we simulated the 
PSF using Quicksirro The simulation shows that 
an r = 30" circle around the position of the pulsar 
contains 88% encircled count fraction for the point 
source. Excluding this circle from the PWN poly- 
gon region (see Figure 3) leaves 3,446 counts in 
MOS1+2, including ~ 60% background contribu- 
tion. Although some regions of the PWN are out 
of the PN Small Window held of view (see Fig- 
ure [5]), we still collected 3,771 counts, including 
~ 70% background contribution. 



Fig. 4. — Smoothed Chandra ACIS image from 
FigureQ]with contours of XMM-Newton MOS1+2 
image overlaid. The contours correspond to the 
brightness levels of 2.9, 3, 3.3, 4, 5, 6, 7, 8, and 9 
counts. 




Fig. 5.— XMM-Newton EPIC PN image is binned 
to a pixel size of 3" and then smoothed with a 
Gaussian kernel of about 9" radius. The image 
also shows the r = 20" circle for the pulsar re- 
gion, 9159 arcsec 2 polygon for the PWN region, 
and r = 60" circle for the background region for 
the spectral analysis of J1357 and its PWN in 0.3- 
10 keV. Unfortunately, the northern part of the 
PWN emission cannot be seen because the source 
is near the chip edge. 



See http://heasarc.gsfc.nasa.gov/docs/xmm/quicksim/quicksim.html 



Table 1: Spectral fits of the field sources in the Chandra and XMM- Newton images 
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Note. — The errors shown represent 90% confidence intervals. 
a PL normalization in units of 10 — 6 photons cm -2 s — 1 keV - x at 1 keV. 
b BB temperature in keV. 

c BB normalization JVr = 10 _3 ^k m /^i0' wnere Kkm i s the source radius in km, and Dio is the distance to the source in units 
oflOkpc. 
d Absorbed flux in the 0.3—8 keV, in units of 10 — 14 erg cm -2 s . 
e Unabsorbed flux in the 0.3-8 keV, in units of 10 -14 erg cm -2 s — 1 . 
Best-fit C-statistic value (source 4) or reduced \ 2 value (other sources). 



2.1.3. Other sources in the field of Chandra and 
XMM-Newton 

Panels a and b of Figure [5] show that Jf 357 is 
the brightest extended X-ray source within the im- 
aged part of the HESS Jf 356 field. Nevertheless, 
we have examined the other sources in the ACIS 
FOV. Three brightest sources on S3 chip south- 
west of the circle in Fi gure^ have Two Mic ron All 
Sky Survey (2MASS; ISkrutskie et ai]|2006h coun- 
terparts. Most likely, they are stars unrelated to 
the TeV source. There are five relatively bright 
sources inside the circle, which do not have 2MASS 
counterparts (Sources 1-5 in Figure Ek). Four of 
these sources are seen in both the ACIS and MOS 
images, while source 2 fell into the chip gap in 
MOS. One more source, lacking a 2MASS coun- 
terpart, is seen in the MOS images (Source 6 in 
Figure [Bb), but it happens to be outside the ACIS 
field of view. We fitted an absorbed PL model to 
the spectrum of each of these sources. When the 
PL model was unacceptable, we fitted the spec- 
trum with the blackbody (BB) model. The spec- 
tral properties of the sources are summarized in 
Table Q] These sources could be CVs or quiescent 
LMXBs, or they could be AGNs seen through the 
Galactic disk. 

Finally, we investigated in detail the relatively 
bright source on ACIS 10 chip, CXOU J135605.9- 
642909 (hereafter CXOU J1356), which is appar- 
ently the same source as 1RXS J135605. 5-642902. 



The source is also seen in the EPIC images, but it 
is near the chip gap in the EPIC MOS1 and MOS2 
images. We extracted 443 total counts within the 
5" radius around the brightest pixel of the source 
from the ACIS image. The spectrum was binned 
to have minimum 25 counts per bin. We did not 
obtain an adequate fit with a single component 
(PL or BB) model (e.g., n K « 2.0 x 10 21 cm- 2 , 
r « 3.4, and xl — 1-8 for 14 degrees of free- 
dom (dof) for the PL model). The spectrum fits 
better the two-component PL+BB model, which 
gives n H « 4.3 x 10 21 cm" 2 , T = 2.9, kT = 0.12 
keV and xl = 1-1 for f2 dof. The absorbed 
flux is 6.4 x f0" 14 erg cm" 2 s" 1 in the 0.3-8 
keV band. We found no significant variability 
in the light curve. We have searched the field 
at other wavelengths to understand the nature of 
the source. The nearest optical/IR source was 
found in the Naval Obser vatory Merged Astro- 
metric Dataset (NOMAD; IZacharias et al.l |2005| ) 
and 2MASS catalogs (B = 13.81, V = 13.66, 
R = 12.69, J = 12.199±0.029, H = 11.779±0.031, 
and K = 11.768 ± 0.030), offset by w 0'.'4 from 
the position of CXOU J1356. The optical-NIR 
colors show that the source might be a G or K 
star. The estimated X-ray to optical flux ratio 
is about 4 x 10 -3 , which is si milar to a typical 
value (~ 10 -3 ) for a G or K star ([Maccacaro et al. 
1988(1 . 

Although many more sources are seen in the 



ACIS and MOS images, the low numbers of counts 
preclude spectral fitting, making it even more diffi- 
cult to establish the nature of these sources. None 
of the investigated sources shows evidence for ex- 
tended emission in X-rays, and there are no com- 
pelling reasons to believe that any of them is re- 
lated to HESS Jf356. 

2.1.4- Multiwavelenths images 

In addition to analyzing the Chandra and 
XMM-Newton images, we have searched the 
archival radio and IR data covering the HESS Jf 356 
region. The 12 micron image (Fig. |BJ:) from In- 
frared Astronomical Satellite (IRAS) shows sev- 
eral bright sources within the circle. Figure [6ji 
shows the 2.4 GHz image from the Parkes survey 
of the southern Galactic planqj- The extended 
radio emission has been prev iously classified as a n 
SNR candidate G309.8-2.6 (JDuncan et al.lll997l) . 
Another radio image from the Pa rkes- MIT-NRAO 
(PMi\_§ iGriffith fc WrightJ Il993h survey at 4.85 
GHz also shows an extended emission southeast 
of HESS J1356 (see Fig. IS:)- The Fermi LAT im- 
age (0.5-100 GeV), shown in Figure 6f, reveals a 
faint GeV emission near the pulsar position. The 
angular extent of the GeV excess is consistent with 
that expected for the unresolved point source. No 
counterparts to J1357 or its X-ray PWN are ap- 
parent in the IR and radio images. Deeper obser- 
vations of this region with higher spatial resolution 
would be desirable to understand the nature of the 
extended radio and IR emission. 



2.2. Spectral and Timing Analysis 

2.2.1. Pulsar Spectrum 

For the spectral analysis of the Chandra obser- 
vation, we first extracted the pulsar's spectrum 
from a small circle (I" radius aperture with 88% 
encircled energy fraction) to reduce the contami- 
nation from the PWN (see Figured]). There are 
371 total counts in this aperture, of which about 
11% come from the PWN and 1% from the back- 
ground. We binned the spectrum with min imum 
of 10 counts per bin and used the C statistic (jCash 
19791 ) in spectral fitting. The spectrum was first 
fitted with a single component model, such as the 



absorbed PL model and the absorbed BB model, 
in 0.3-8 keV. The PL fit gave n H . 2 i=n H /{l0 21 
cur 2 ) w 1.4, and T w 2.4 (C = 52 for 34 bins, 
31 dof ), while the BB was not statistically ac- 
ceptable (C = 173 for 34 bins, 31 dof ). Since 
the PL fit quality was not good, we fit the same 
spectrum with the absorbed PL+BB model. The 
much better fit (C = 28 for 34 bins for 29 dof 
) yields n H ,2i « 4.7, T « 1.7, kT « 0.14 keV 
(T f» 1.6 x 10 6 K), and the projected emitting area 
A as 13 d 2 . 5 km 2 (d 2 .g = d/2.5 kpc). This temper- 
ature and area correspond to the apparent radius 
R = sJApK ~ 2d 2 .5 km and bolometric luminos- 
4ArT 4 - 2 x 10 32 d2. „ erg s" 1 . The 
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encircled energy fraction), while the correspond- 



unabs 
psr 



= 32 



+ 190 
-21 



x 10 
The un 



14 



10 32 d 2 5 erg 



ing unabsorbed flux is F } 

erg cm^ 2 s _1 in the same energy band. 

absorbed luminosity is L psr = 2.4 

s _1 . The luminosity of the non-thermal (PL) 

component is Lpl = 5.2t 2 ;g x 10 31 erg s _1 (~ 

22% of the total unabsorbed luminosity), while 

the luminosity of the thermal (BB) component is 

L B b = l-9^i 4 4 x 10 32 erg s _1 . The uncertainties of 

the fits are reported in Tabled and the confidence 

contours are shown in Figures [7] and HI 

Instead of BB, one can use a neutron star 
atmosphere (NSA) model for the thermal com- 
ponent, which gives a lower (effective) temper- 
ature and a large r size of the emitting region 



(jPavlov et al.lll995l ). Since the spectrum emerging 



from a magnetized atmosphere is anisotropic, one 
should know the geometry of the surface magnetic 
field (e.g., the angles between the rotation axis, 
magnetic axis and line of sight for a dipole field) 
to fit the spectrum. For a crude estimate, however, 
one can use the NSA models in XSPEC, which are 
calculated for emergent flux assuming the mag- 
netic field perpendicular to the surface. We made 
such an estimate using the hydrogen NSA mod- 
els for the NS mass M = 1.4M , radius R = 10 
km, and magnetic field B = 1 X 10 13 G. If the 
distance is fixed at d = 2.5 kpc, the fit (C = 29 
for 30 dof) yields ra H ,2i = 4.6+^, T = 1.45l^°, 
T cff = 0.96±g;gt MK (T£ = 0.74l°;° 3 MK for the 
gravitationally redshifted temperature Jj. It shows 



7 see http://www.atnf.csiro.aU/research/surveys/2.4Gh_Southern/datMhtSJItilar fit 
8 see http://www. parkes. atnf.csiro.au/observing/databases/pnin/pnin.'rrtrnty 



obtained by iLemoine-Goumard et all 




Fig. 6. — X-ray, IR, and radio images of the J1357 and HESS J1356 field. The cross shows the centroid 
position of HESS J f 356, while the circle (15' radius) indicates the intrinsic (rms) Gaussian width of the source, 
according to lRenaud et al.l ( 20081 ). The white arrows in panels (a) and (b) indicate CXOU J135605. 9-642909 
and 6 other sources (see text for details), (a) Chandra ACIS image in 0.5-8 keV. (b) XMM MOS1+2 image 
in the 0.3-10 keV. (c) 12 micron IRAS image, (d) 2.4 GHz Parkes image. The small circle located south of 
the cross point indicates the position of SNR candidate G309.8-2.6. (e) 4.85 GHz PMN image, (f) Fermi 
LAT image. The diamond in panels c-f shows the position of J1357. 
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Fig. 7. — Top: Confidence contours (68% and 
90%) in the T - n H plane for the PL+BB fit to 
the pulsar's spectrum from the Chandra observa- 
tion (1" extraction radius) and Chandra + XMM- 
Newton observations (20" extraction radius). Bot- 
tom: Confidence contours (68% and 90%) for the 
PL component of the PL + BB fit to pulsar's spec- 
trum from the Chandra observation (1" extrac- 
tion radius) and Chandra + XMM-Newton obser- 
vations (20" extraction radius). The PL normal- 
ization is in the units of 10~ 5 photons cm~ 2 s _1 
keV -1 at 1 keV. The dashed curves are the lines 
of constant unabsorbed flux for the PL component 
in the units of 10 -14 ergs cm -2 s _1 in 0.5-8 keV. 



that the thermal component can be emitted from 
a substantial part of the NS surface. 

To compare the Chandra ACIS spectrum with 
those obtained from the XMM-Newton EPIC de- 
tector, we extracted the pulsar spectrum from the 
20" radius aperture for MOS and PN and the 



Fig. 8. — Confidence contours (68% and 90%) for 
the BB parameters of the PL + BB fit to the 
pulsar's spectrum from the Chandra (1" extrac- 
tion radius) observation and Chandra + XMM- 
Newton observations (20" extraction radius) . The 
BB normalization is the projected emitting area 
(km 2 ). The dashed lines show the constant bolo- 
metric luminosity in units of 10 32 ergs s _1 , assum- 
ing d = 2.5 kpc. 

background spectrum from the 100" radius for 
MOS and 60" radius for PN (see Figures [3] and 
[5|). The spectrum is heavily contaminated by the 
compact PWN (about 42% contribution). The 
spectra were binned with minimum 50 counts per 
bin for each detector. Since the absorbed PL fit 
was not satisfactory (x 2 = 1.94 for 82 dof), we 
fitted the absorbed PL+BB model which yields 
n H ,2i ~ 2.2, T « 1.35, kT w 0.18 keV, and the 
projected emitting area A « 2d 2 . 5 km 2 . The 
BB parameters correspond to the apparent ra- 
dius R ~ 0.8d 2 .5 km and L hoi = 8.6 x 10 31 d| 5 
erg s _1 . Although the best-fit kT and A from 
the XMM-Newton observations are different from 
those observed from the Chandra data, these pa- 
rameters are consistent within their uncertainties. 
The observed flux (corrected for the 80% encircled 
counts fraction) is F^ aal+PWN = 2.6 x 10~ 13 erg 
cm _2 s _1 , and the corresponding unabsorbed flux 

i« ^pTlsar+PWN = 3-6 X 10 



" 13 erg cm 2 s l in 0.3- 
10 keV. Additional details are provided in Tablc[2] 
The photon index and temperature of the 
PL+BB model inferred from the XMM-Newton 
and Chandra observations are consistent within 
their uncertainties. On the other hand, the ab- 



Table 2: Absorbed BB+PL fits to the pulsar (or pulsar + compact PWN) spectrum from the Chandra and 
XMM-Newton observations. 



Model 3 



"■H,21 



AP 



kT c 



A d 



7 +3. b 
-2.8 
22+ i.o 



,+U.bb 
0.63 
i qcr + 0.16 

-0.50 
-0.69 

,+0.16 
0.17 



TTTT- 
5.8 
22 ±5 



+0.Ub 

0.04 
+0.04 

0.03 
+ 0.13 

0.09 
+0.04 

0.03 



+4b4 
12 
n+6 
z -l 
+22 
0.4 
+2.8 
0.6 



- abs 



qrg — 

-0.8 

17.8±0.8 

16 ±2 
17.3 ±0.7 



? — iT 7 " 

unabs 



+iyo 

21 
+8 

-4 
+31 



(C or X j)/dof s 



Chandra (1") 

XMM-Newton (20") 

Chandra (20") 

Combined (20") 



4.7^ 



2.4 
1.9 



0.7 
+3.1 

1.9 
+0.8 

0.6 



1.72"* 



1.32Z 
1.33" 1 



0.8 



18 



+ 18 
-12 



21 ±5 



0.14 
0.18 
0.22 
0.19 



13 

2 

0.5 

0.9 



6.7^ 



32 
26 
22 



28.2/29 
1.27/80 
1.03/82 
1.18/167 



Note. — The errors shown represent 90% confidence intervals. 
Radius of extraction region is given in parentheses. 



b Spectral flux in units of 10 
c BB temperature in keV. 



photons cm 2 s 1 keV 1 at 1 keV. 



d Projected emitting area for BB model in units of km 2 at the distance of 2.5 kpc 



Absorbed flux in the 0.3-8 keV, in units of 10 -14 



erg cm 



Unabsorbed flux in the 0.3—8 keV, in units of 10 14 erg cm 2 s 1 . 
g Best-fit C-statistic value for Chandra observation (1" radius aperture) or reduced \ 2 value (other cases). 



sorbed flux estimated from the XMM-Newton 
data is a factor of three larger than that esti- 
mated from the Chandra data due to the larger 
extraction region (20" radius), chosen for XMM- 
Newton data, which contains a significant PWN 
contribution. Therefore, we extracted the spectra 
using the same source region (20" radius) for the 
Chandra observation to compare with the spectral 
results obtained from the XMM-Newton observa- 
tion. The region includes 873 source counts, which 
is 83% of the total counts. Now both fits are con- 
sistent with each other (see Table. [2]). Finally, the 
Chandra (20" radius aperture) and XMM-Newton 
data were fitted jointly with the PL+BB model, 
which still fits satisfactorily while the uncertainties 
of the best-fit parameters are significantly reduced 
compared to the individual fits (see Table [2] and 
Figured]). 








1 2 5 

Energy (keV) 



Fig. 9.— ACIS-I (black), EPIC MOS1 (red), 
MOS2 (green), and PN (blue) spectra of the pul- 
sar region (20" radius) fitted with the absorbed 
PL (dashed line) -I- BB (dash-dotted line) model. 



2.2.2. PWN Spectrum 

We extracted the ACIS spectra from the PWNI 
and PWN2 regions (see Fig. [J). For the PWNl's 
spectrum, we excluded the pulsar region (1"5 ra- 
dius aperture) and used the absorbed PL model in 
the energy range of 0.3-8 keV. The spectrum was 
binned with minimum of 10 counts per bin. The 
fit (xl = 1-3 for 24 dof) gives n H ,2i = 3.2±?;|, T = 
(6.0±0.9)'xl0- 14 ergcm- 2 



h0.37 



i oi+u.oi j pabs 

1 - O1 -0.33' dIlu r pwnl 

s -1 . The spectrum of PWN2 was also binned with 
minimum of 10 counts per bin and fitted with the 
absorbed PL model in the same energy range. We 



used the C-statistic because of small number of 
counts in the region. The PL fit to the spectrum 
of the PWN2 gives somewhat larger T = 1.8to'g; 
however, the uncertainties are large as one can 
see from the confidence contours shown in Fig- 
ure [TUl We also extracted and fitted the spectrum 
from both regions combined, since the difference 
in slopes of both spectra is not statistically signifi- 
cant. The fit gives u H ,2i = 3.7±f%, T = 1.37±0.22, 
and the absorbed flux is F^ s nl+2 = (7± 1) x 10" 14 
erg cm _2 s _1 in the 0.3-8 keV band, while the 



10 



unabsorbed flux is F™^_, 



*t? 



10 



-14 



erg 



cm" s~ in the same energy band. It corresponds 
to the luminosity -LpwNi+2 ~ 7 x 10 31 d 2 5 erg s" 1 , 
i.e., about 130% of the pulsar's non-thermal lumi- 
nosity estimated from the same observation. 




Fig. 10. — Confidence contours (68% and 90%) in 
the «h — T plane for the PWN spectrum measured 
with XMM-Newton (blue) and the PWN1 {black), 
PWN2 (red), and PWNl±2(green) spectra mea- 
sured with Chandra. 

For the PWN region on the II chip (see Fig.[lJ), 
we binned the spectrum with minimum 25 counts 
per bin and fitted with the absorbed PL model 
in the 0.3-8 keV energy range. We fixed the hy- 
drogen column density to the value tih,21 = 4.7, 
found from the PL+BB model for pulsar obtained 
from the Chandra data (see Table 2). The fit 
(xl = 0.66 for 10 dof) gives T = 2.5±^;| and 
7V r = 8±l x 10" 6 photons cm" 2 s" 1 keV" 1 at 1 
keV. The absorbed flux in the same energy range 
is 1.3±g;? x 10" 14 erg cm 
sorbed flux is 41}° x 10" 14 
spectrum of the PWN region on the 12 chip was 
binned with minimum of 50 counts per bin and 
fitted with a single absorbed PL model in the en- 
ergy of 0.3-8 keV. We fixed the hydrogen column 
density of n.n,2i = 4.7. The fit was acceptable 
(xl = 0.7 for 18 dof) with T = 1.7 ± 0.4 and 
A/r = (1.5 ±0.5) x 10" 5 photons cm" 2 s" 1 keV" 1 
at 1 keV. The absorbed flux in the same energy 
range is (7 ± 3) x 10" 14 erg cm" 2 s , while the 
unabsorbed flux is (9 ± 2) x 10" 14 erg cm" 2 s" 1 . 

Figure [TT1 shows the MOS1 and MOS2 spectra 
of the region shown in Figure [3l binned with mini- 



is 1.3^07 x 10 14 erg cm 2 s 1 , while the unab 



erg cm 2 s 1 . 



The 



mum 50 counts per bin. The region is much larger 
than that of PWN1+PWN2 of Chandra (see Fig- 
ures UJ and [3]). The absorbed PL fit is good (xl 
= 0.94 for 61 dof). The best-fit hydrogen column 
density and photon index are nn,2i — 3.9j^;f and 



1.7 



j^abs 
pwn 



-0.3 

-0.2' 



respectively. The absorbed flux is 



(5.0 ±0.6) x 10 



-13 — 9 - 

erg cm s 



in 0.3- 



8 keV. The unabsorbed flux F, 



unabs 



7.5 



+ 1.4 



10 13 erg cm 2 s 



pwn — ' -"-0.9 ~ 

The results of the fits, sum- 



marized in Table |31 suggest that the PWN spec- 
trum softens with the distance from the pulsar. 
The combined compact±extended PWN flux is 
-fpwn bs ~ 9 x 10" 13 erg cm" 2 s" 1 which corre- 
sponds to the the luminosity L pwn 



6.7xl0 32 d| 5 



erg s 



2.8L L 




rrs 
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Fig. 11.— EPIC MOS1 (black) and MOS2 (red) 
spectra of the PWN shown in Figure [3) 

2.2.3. Timing Analysis 

We searched for pulsations in the PN data (time 
resolution of 6 ms) extracted from the r = 20" cir- 
cle centered at the pulsar position, which contains 
939 counts, including 15% from the background, in 
the 0.3 - 1.1 keV band, where the BB component is 
dominant (see Figure [9]). From the same region in 
the 1.1 - 10 keV band, we extracted 1473 counts, 
including 21% from the background. The pho- 
ton arrival times were corrected to the solar sys- 
tem barycenter using the SAS barycen tool. The 
PN observation started at the epoch 55,037.63605 
MJD and continued for the time T span = 55.5 ks. 
We used the r adio t iming ephemeris reported by 
' 20041) (/ = 6.020167 772 6(4) s" 1 



Camilo et al 



and / = -1.305 395(4) x 10" 11 s" 2 for the epoch 
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Table 3: Spectral properties of the PWN from the Chandra and XMM-Newton observations 



Region 


1H.21 


r 


Af a 


-**abs 


-^unabs 


(C 


or x^)/dof d 


Chandra 
















PWN1 


n q+2.3 
au -XA 


-, O1+0.37 
- 1 — 31 -0.33 


8 5+ 4 ' 7 
0.0_2 g 


6.0 ±0.9 


7±1 




1.3/24 


PWN2 


6-o^; 9 5 


1 s +1 -° 


q W+9.5 

°-'-3.7 


1 O + 0.5 

1 -°-0.4 


n 9+4.2 
z - z -0.8 




1.04/5 


PWN 1+2 


o 7 +2.0 
"'•'-1.3 


1.37 ±0.22 


12±S 


7±1 


9+ 2 




0.90/31 


PWN 11 


4.7(fixed) 


2 5+ 1 ' 8 


8±S 


i ,+0.9 
J -°-0.7 


4 +io 




0.66/10 


PWN 12 


4.7(fixed) 


1.7 ±0.4 


15 ±5 


7±3 


9±2 




0.70/18 


XMM-Newton 
















PWN 


3-»ii:i 


1 7 4+0.27 


123±« 


50 ±6 


75t^ 4 




0.94/61 



Note. — The errors shown represent 90% confidence intervals. 
a Spectral flux in units of 10 -6 photons cm -2 s —1 keV~ 1 at 1 keV. 
b Absorbed flux in the 0.3-8 keV, in units of 10~ 14 erg cm -2 s" 1 . 
c Unabsorbed flux in the 0.3-8 keV, in units of 10 -14 erg cm -2 s — 1 . 
" C-statistic value for PWN2, x 2 values for other regions. 



of 52,921.0 MJD) to calculate the predicted fre- 
quency at our epoch. There was a glitch at the 
epoch of 52,021±16 MJD (ler uncertainties), but 
it was before the ephemeris was measured. We 
assumed that there were no glitches after that. 
The predicted frequency for our observation is 

/predict - 6.017 780 858(7) Hz. 

We applied the Z\ test (jBuccheri et al.lll983 ; 



Pavlov et al.l Il999l) to search for the pulsed sig 



nal in the band 6.017 780 ± 0.000 010 Hz near the 
predicted frequency, with a step of 0.1 ^/Hz w 
0.006T S pa n . We found the maximum Z^ ^ = 40.5 
at / = 6.017 7801 Hz ± 1.6 ^Hz, i.e.,' the pulsa- 
tions are detected with a 6er significance in the 
0.3-1.1 keV energy rangq^j. The maximum value 
of Z\ in the 1.1-10 keV energy range is only 7.9 at 
/ = 6. 017 78 19 Hz ± 3.5 /xHz, which corresponds 
to 2.3cr significance. 

Using the two epochs (55,037.63605 MJD from 
our observation and 52,921.0 MJD from the ra- 
dio) and two frequency values, we can estimate 
the frequency derivative, / = —1.3056(7) x 10~ n 
s~ 2 , which is consistent with / determined from 
the radio. Figure [T2l shows the pulse profiles with 
10 phase bins, with the estimated pulsed frac- 
tion (the ratio of the number of counts above 



■'The la frequency uncertainty is calculated as Sf = 
3 1 / 2 7r- 1 T sp L(Z 2 imax )" 1/2 . This equ ation can be derive d 
using the method outlined in §2.4 of iBretthorst I (1988). 
The numerical factor, V3/n ~ 0.55, was confirmed in 
Monte Carlo simulations by A. Kaurov (2011, priv. comm.). 



the minimum level to the total number of counts) 
p = 36% ± 5% in 0.3-1.1 keV, 18% ±4% in 1.1-10 
keV, and 14% ± 3% in 0.3-10 keV. We es timated 
the e rror of pulsed fraction as Sp = y/Z/N (|Linslev 
1975T ). After correcting for the background, we ob- 
tain the intrinsic pulsed fraction p; nt = 42% ± 5%, 
23% ± 4%, and 17% ± 3%, respectively. Thus, we 
conclude that the X-ray emission from J1357 is 
undoubtedly pulsed, althoug h the pulsed fraction 
is not as high as reported by IZavlinl (|2007l) . Inter- 
estingly, the pulsed fraction is higher at lower pho- 
ton energies, where the thermal component domi- 
nates. Such strong pulsations of the thermal emis- 
sion are likely due to the anisotropy of atmospheric 
radiation in the strong magnetic field of this pul- 
sar. 

3. Discussion 

Our Chandra ACIS and XMM-Newton EPIC 
observations have shown that the PWN of J1357 is 
much more extended and luminous than it could 
be seen from the earlier observations. We have 
also found that the thermal component of pulsar 
emission is strongly pulsed, while the pulsation of 
the nonthermal component, which dominates at 
higher energies, are much weaker. Below we dis- 
cuss the physical implications of our findings, com- 
pare the J1357 PSR/PWN properties to those of 
other similar pulsars observed in X-rays, and in- 
vestigate the relation between J1357, HESS J1356, 
and putative host SNR G309.8-2.6. 
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Fig. 12. — Pulse profiles of J1357 extracted from 
the EPIC PN data in three energy bands. 

3.1. Pulsar 

The pulsar spectrum cannot be satisfactory fit- 
ted by simple one-component models, such as ab- 
sorbed PL or BB (see §2.2.1). The two-component 
BB+PL model fits the spectrum much better, al- 
though some residuals are still present in the com- 
bined fits. The best-fit photon index for the PL 
component, F ~ 1.7±0.6, is marginally larger than 
T ~ 1.3±0.3 of the compact PWN (region PWN1). 
The non-thermal luminosity is L non th ~ 5 x 10 31 
erg s _1 in the energy range 0.3-8 keV. The PL 
component may partly come from an unresolved 
PWN. This could explain the weak pulsations 
above 1 keV (see §2.2.3), in contrast to, for in- 
stance, the Vela pulsar (which is well resolved 
from the surrounding PWN), whose pulsed frac- 
tion grows with en ergy reaching 62% above 1.8 
keV (|Sanwal et alj 120021) . Our timing analysis 
limits the luminosity of the unresolved PWN to 
< 77% of the point source luminosity in 1.1-10 
keV. 



Unlike the hard non-thermal component, the 
soft thermal component should not be strongly 
contaminated by the PWN emission since it 
clearly dominates at < 1 keV. Therefore, the best- 
fit BB temperature, T w 1.9 - 2.7 MK, radius, 
R ~ 0.3 — 1.1 km, and bolometric luminosity, 
iboi ~ 2x 10 31 erg s _1 , are reasonably accu- 
rately determined (see Table [2|) . The values are 
similar to those of other Vela-like pulsars. Al- 
though the temperature is a factor of two higher 
than expected from the standard NS cooling, and 
the emitting area radius is a factor of ten smaller 
then the classical NS radius, the discrepancies 
can be attributed to the use of the simplistic 
BB model. Indeed, the f it with the hydr ogen 



atmosphere models (e.g.. iPavlov et al.l 119951 ). de- 



scribed in Section 2.2.1, gives a lower temperature, 
T™ = 0.70-0.78 MK at R°° = 13 km, at the same 
bolometric luminosity. The measured bolometric 
luminosity is a factor of ~ 3 lower than predicted 
by the theoretical cooling curve for low-mass NSs 
(Mns = 1.3M©) without superfluidity; it is con- 
sistent with the cooling curves for heavier NSs 
'e.g., Mns = 1.5-1.6 M&) for som e superfluidity 



models (JYakovlev fc Pethickl 120041 ) . 

Very surprising is the high (w 40%) pulsed frac- 
tion associated with the thermal component. In 
general, one should not expect such strong pul- 
sations from a non-uniformly heated NS surface 
with dipolar magnetic field for a blackbody-like 
angular distribution. However, there have been 
recent reports of even stronger thermal pulsa- 
tions in strongly magnetized young pulsars, such 
as PSR J1119-6127 (74% ± 14%; Gonzalez et al. 
2005) and PSR J1718-3718 (52% ± 13%; Zhu et 
al. 2011). This suggests that the strong pulsa- 
tions are caused by the combined effect of the 
anisotropic temperature distribution and strong 
anisot ropy of emitted r adiation in high magnetic 
fields ()Pavlov et al.lll994l ) . An alternative possibil- 
ity is that part of emission interpreted as thermal 
from our BB+PL and NSA+PL fits is in fact an 
additional, softer non-thermal component, which 
is masked at higher energies by the harder, unre- 
solved PWN emission. This unaccounted compo- 
nent might be the reason for the remaining resid- 
uals of the PL+BB and NSA+PL fits. However, 
higher quality data are needed to test this hypoth- 
esis. 
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3.2. PWN 

At the plausible distance of 2.5 kpc, the 
unabsorbed X-ray luminosity of the compact 



PWN, L 



pwnl 



5 x 



10 31 ergs s" 1 in the 0.3- 



keV band, corresponds to the X-ray efficiency, 
77 pwnl = L pwn i/E ~ 1.7 x 10 -5 , lower than those 
of most Vela-like pulsars (see Kargaltsev et al. 
2007; KP08). This suggests that either the dis- 
tance is underestimated or, in addition to E and 
t, ?7 P wni strongly depends on other factors (e.g., 
the pulsar's speed and the angle between the spin 
and magnetic axes). Note, however, that the effi- 
ciency of the entire detected PWN is higher by a 
factor of 15, ry pwn ~ 2.6 x 10~ 4 . 

The spectral slope of the compact PWN, 
Tpwni = 1.3 ± 0.3, is similar to those of PWNe 
around Vela-like puls ars (listed in Table 2 of 
iKargaltsev et aJjl2007h . The PWN spectrum ap- 
parently softens with the distance from the pul- 
sar, becoming V sa 1.7— 1.8 for the narrow feature 
(PWN2) and the large-scale PWN, which suggests 
synchrotron cooling of the outflow. The degree of 
spectral softening is similar to that seen for the 
extended tail of PSR J1509-5850, where the pho- 
ton index changes from V — 1.8 ± 0.3 for PWN 
in the vicinit y of the pulsar to 2.4 ± 0.4 for the 
extended tail (JKargaltsev et al.ll2008l ). 

Our X-ray images reveal the morphology of the 
PWN on different angular scales. The observed 
morphology does not easily fit in either bow-shock 
or torus-jet category (see KP08). The bright com- 
pact PWN1 extends up to ~ 10" from the pulsar 
and has rather amorphous, somewhat asymmet- 
ric morphology. The diffuse emission is notice- 
ably brighter northeast of the pulsar compared 
to the opposite side, although some emission is 
still clearly discernible within the PWN1 region 
southwest of the pulsar. We find no evidence of 
structures that could be interpreted as a torus as- 
sociated with the termination shock (TS) in the 
pulsar wind, which may suggest that the TS oc- 
curs at an angular distance of several arcseconds 
or less, corresponding to r s < 4x 10 16 d2.5 cm (ver- 
sus 1 x 10 17 cm in the Vela PWN; Hclfand et al. 
2001). If the pulsar is moving subsonically inside 
the hot SNR medium, we can use the above esti- 
mate to obtain a lower limit on the ambient pres- 
sure at p amb ~ E^ircr^y 1 > 5 x 10 _9 (r s /4 x 
10 16 cm) -2 dyn cm -2 , which is close to the maxi- 



mum pressure one can expect in a 10 kyr ol d SNR 
(JKargaltsev et al.ll2009t iBamba et al.ll2010h . The 
high pressure required, and the asymmetry of the 
compact PWN, elongated approximately along the 
northeast-southwest direction, might be indicative 
of the TS being also asymmetric, i.e., being far- 
ther away from the pulsar on the northeast side 
than on the southwest side. Such situation may 
occur if the pulsar is moving fast through the 
surrounding medium in the southwest direction, 
and the ram pressure of the medium compresses 
the PWN in front of the pulsar. If the pulsar's 
speed exceeds the sound speed in the medium, 
c s = (5fcT/3A*TO H ) 1/2 = 12 n~^ 2 Tl /2 km s~\ a 
bowshock is formed. The appearance of a bow- 
shock PWN depends on the Mach number, uni- 
formity of the medium, and intrinsic anisotropy of 
the pulsar wind. X-ray PWNe around superson- 
ically moving pulsars often exhibit a bright "bul- 
let" in the vicinity of the pulsar and a much more 
extended faint tail in the direction opposite to 
that of the pulsar motion; however, more complex 
structures have been seen in some cases (e.g., the 
Guitar and Eel PWNe associate d with B2224+65 
and J1 826-1256, resp ectively; Ijohnson fc Wang 
(l2010h ; lRoberts et alJ (l2007ll l. 

The adjacent PWN2 region shown in Figure 1 
encompasses the fainter narrow feature attached 
to the compact PWN. Interestingly, this narrow 
feature bends sharply at about 23" from the pul- 
sar. Although there is no doubt that the feature 
is part of the PWN, the origin of the feature is 
not clear. For instance, it could be a pulsar jet. 
Indeed, some of pulsar jets are known to show ex- 
tre me bending (e. g ., the outer jet of the Vela pul- 
sar; IPavlov et al.1 ( 2003 )). The apparent lack of 
a similarly looking counter-jet should not be sur- 
prising since the two jets can differ significantly 
both in shape and surface brightness due to the 
Doppler boost and/ or proper motion effects (e.g., 



Pavlov et al 



(J2003I )). On the other hand, the ra- 
tio of the compact PWN luminosity to that of the 
putative jet, £ pwn i/L pwn 2 ~ 5, is noticeably larger 
than that in the Vela and Crab PWNe but com- 
par able to that of the P WN around PSR B1706- 
44 (JRomani et al.l 120051) . If this interpretation is 



taken at the face value, one could expect the pulsar 
to be moving northeast (in general, jets tend to be 
co-aligned with the pulsar's di r ection of motion-cf. 
the Vela pulsar; IPavlov et all (|2003h ). 
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Alternatively, the narrow feature could be a 
part (e.g., an inner channel ) of the pulsar tail (see, 
e.g., iKargaltsev fc Pavlovl (|2008l ) and references 
therein). In this case the pulsar would be mov- 
ing in the opposite (southwest) direction. This 
interpretation is in line with the above-described 
asymmetry of the compact PWN whose brighter 
part could, in this case, indeed be the "bullet" 
associated with the termination shock deformed 
by the ram pressure of the oncoming ISM. How- 
ever, the sharp bending and the lack of a fainter, 
much more extended tail (cf. J1509-5850 tail; 
Kargaltsev &; Pavlovl (|2008l )). and the emission 
west-southwest of the compact PWN, are at odds 
with this interpretation 1 H 

Finally, the narrow feature could be a part of 
the bow associated with the forward shock (see 
e.g., the Geminga PWN; Pavlov et al. 2006, 2010) 
with the pulsar moving southeast. In this case 
the observed large-scale emission can also be as- 
sociated with the asymmetric bows hock similar to 



that seen around PSR J1826-1256 (Ro berts et al. 



20071 ). Somewhat puzzling may be the lack of any 
X-ray emission within the bow interior; however, 
this is also the case in the Guitar and Eel bow- 
shocks. To conclude, with the data in hand we 
cannot univocally interpret the morphology of the 
X-ray PWN. Measuring the pulsar's proper mo- 
tion should provide decisive information allowing 
one to discriminate between the above interpreta- 
tion s. Assuming a typi cal pulsar speed of 400 km 



s " (jHobbs et al.1120051) , this should be possible to 
accomplish in several years with VLBI. 

3.3. The Nature of the VHE source 

Given the lack of other promising counterparts 
and the location of J1357 within the HESS J1356 
extent, it seems very plausible that the two 
should be related in some way. Yet the nature 
of HESS J1356 has not been firmly established. 
The TeV emission could be attributed to either 
the host SNR or to the relic PWN of J1357. To 
date, several well-known SNRs have been detected 



in the TeV band (e.g.. lAharonian et al.ll2004l ). In 
the cases of resolved SNRs, the TeV emission was 



The diffuse emission seen in the PMN 4.85 GHz image 
north of the pulsar (panel e in Fig. 6) could, in principle, be 
the radio counterpart of the extended pulsar tail. Deeper 
radio observations are required to test this hypothesis. 



associated with the SNR shell, and the sizes of 
the SNRs in TeV images were similar to those of 
the non-thermal radio shells indicating that both 
emission components are powered by particles ac- 
celerated in the forward shock. Although no shell 
is seen in the radio images of HESS J1356, the 
radio emission is mainly seen near the center of 
HESS J1356 (see Fig. \$. Thus, in this case the 
TeV emission does not appear to be associated 
with the nonthermal shell, and the shell itself is 
not seen in the radio. 

An alternative, more plausible, interpreta- 
tion of the TeV and radio emission could be a 
relic PWN. Many extended TeV sources neigh- 
bor young Vela-like pulsars, sometimes offset up 
to 10'-20' from the center of TeV emission (see 
KP10 for recent review). Furthermore, recent X- 
ray observations of the Vela pulsar region (Vela 
X; Mori et al. 2008), PSR J1826-1334 (Gaensler 
et al. 2003; Pavlov et al. 2008) a nd PSR J1809- 
1917 (JKargaltsev fc Pavlovl 120071 ) have provided 
convincing evidence that the TeV sources are con- 
nected to the parent pulsars through faint asym- 
metric X-ray nebulaa 12 ! (see Table 0]). Indeed, 
the Chandra observations have demonstrated that 
in the above three examples the PWNe consist 
of compact (0.1-0.8 pc) bright cores and fainter 
asymmetric, more extended (2-8 pc) components. 
The offsets and the asymmetries of the X-ray 
PWNe could be created by the reverse SNR shock 
that had propagated through the nonhomoge- 
neous SNR interior and reached one side of the 



PWN sooner than the other side (Blondin et al 
120011 ). This scenario could also account for the 
similarly asymmetric, offset TeV emission (e.g., 
de Jager & Djannati-Ata'i 2009). However, the 
physical origin of the TeV emission still remains 
under debate. It can be produced via the In- 
verse Compton Scattering (ICS) of the relativistic 
pulsar-wind electrons on photons from the om- 
nipresent cosmic microwave background (CMB), 
galactic IR background, and IR photons from local 
star-forming regions and warm dust clouds. Al- 
ternatively, the TeV photons can be produced as a 
result of 7T° — > 7+7 decay, with tt° being produced 
when the relativistic protons from the pulsar wind 



12 There remains a large fraction of relic PWN candidates 
where only bright compact PWN is seen X-rays while the 
extended component is not seen perhaps due to is faintness 
and/or rapid cooling of the outflow. 
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Table 4: Properties of candidate TeV plerions and their parent pulsars. 



HESS ID 


A' 


TcV size 


Offsct b PSR 


£ 36 c 


: T 


Dist. 


L^/E 


Lx/L~f 


R • d 






pc 








kyr 


kpc 








J1809-193 


0.3 


20 


8' 


J1809-1917 


1.8 


50 


3.5 


1% 


2% 


8 


VelaX 


0.75 


5 


30' 


J0835-4510 


6.9 


11 


0.3 


0.002% 


90% 


- 


J1825-137 


0.48 


70 


10' 


J1826-1334 


2.8 


21 


3.5 


4% 


0.3% 


7.5 


J1356-645 


0.79 


40 


8' 


J1357-6429 


3.1 


7.3 


2.5 


0.2% 


10% 


19 



Unabsorbcd 7-ray flux (1-10 TeV) in units of 10 



-11 



ergs s cm 



b Offset of the TeV source center from the pulsar. 

c Pulsar spin-down power, E, in units of 10 36 ergs s — 1 . 

d Average surface brightness measured from Chandra ACIS images, in counts ks" 



interact with the ambient matter (e.g., Horns et 
al. 2007). In the case of HESS J1356, the offset 
is well within the range observed in other relic 
PWN candidates (see KP10). The ratio of the 
TeV luminosity to the pulsar's spin-down power, 
LTev/E ~ 0.002, and the ratio of X-ray to 7-ray 
(1-10 TeV) PWN luminosities, L x /L TcV ~ 0.1, 
are similar to those of other relic plerions. 

Alternatively, the relic PWN could be left be- 
hind the fast moving pulsar if the latter is moving 
in the northeast direction (see above). At a typi- 
cal speed of 400 km s _1 , the pulsar would move by 
5.'6 in 10 kyrs, which is consistent with the offset 
between the pulsar and the center of HESS J1356. 
So far, there have been no reports of VHE emission 
from pulsar tails. It would be interesting to estab- 
lish the presence of the TeV emission from an ex- 
tended pulsar tail because this would help to break 
the degeneracy in interpreting the nature of the 
VHE emission from the crushed PWNe, for which 
both leptonic (Inverse Compton) and hadronic (ir° 
decay) TeV emission mechanisms are currently be- 
ing debated. If the n° decay is the dominant pro- 
cess in the crushed PWNe, we would expect pul- 
sar tails to be significantly fainter in the TeV be- 
cause high-speed pulsars move in low density me- 
dia outside their host SNRs. On the other hand, 
in the leptonic scenario, crushed plerions and pul- 
sar tails should have comparable TeV luminosities. 
Regardless of the nature of the relic PWN and the 
mechanism responsible for the VHE emission, one 
can expect that the observed extended radio emis- 
sion within HESS J1356 may also be produced by 
the pulsar wind if the electron SED extends to suf- 
ficiently low energies. Better quality radio images, 
more complete multiwavelength spectrum, and ro- 



bust multizone pulsar wind models are required to 
assess the origin of the radio emission. 
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